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Abstract 

We introduce in this paper two dimensional lattice models whose continuum limit belongs 
to the N — 2 series. The first kind of model is integrable and obtained through a geometrical 
reformulation, generalizing results known in the k — 1 case, of the Tk vertex models (based on the 
quantum algebra U q sl(2) and representation of spin j = k/2). We demonstrate in particular that 
at the N — 2 point, the free energy of the Tk vertex model can be obtained exactly by counting 
arguments, without any Bethe ansatz computation, and we exhibit lattice operators that reproduce 
the chiral ring. The second class of models is more adequately described in the language of twisted 
N = 2 supersymmetry, and consists of an infinite series of multicritical polymer 

points, which should lead to experimental realizations. The presence of N = 2 in that case 
is traced back to the Parisi Sourlas supersymmetry of the lagrangians usually used to replace 
n — > limits. Boundary conditions as well as fermionic and bosonic variables are geometrically 
interpreted. Moreover it turns out that the exponents v = (k+ 2)/2(fc + 1) for these multicritical 
polymer points coincide with the old phenomenological formulas of Flory. We therefore confirm 
that these formulas are exact in two dimensions, and suggest that their unexpected validity is due 
to non rcnormalization theorems for the N = 2 underlying theories. We also discuss the status of 
the much discussed theta point for polymers in the light of N = 2 rcnormalization group flows. 
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1 Introduction 



N = 2 theories are one of the most interesting example of conformal theories in two dimensions. 
A large subset of them possess in particular a very useful Landau Ginzburg description p], 0, Bj , 
whose efficiency is ensured by non renormalization theorems Q ^. The Landau Ginzburg description 
makes also transparent the connection with singularity theory 1 1, , and explains nicely the origin 
of ADE type classification. 

Unfortunately these beautiful structures have not so far been much observed. It is known in 
principle how some special points of the Tk vertex models the XY (Gaussian) model or the 
Ashkin Teller §, @, § model possess N = 2 supersymmetry. However it has not been clear so far 
what is special about these points, and therefore how to tune the parameters to observe them. 
Also one would like the Landau Ginzburg description, and the chiral ring, to play some more 
illuminating role in these identifications, and make sense physically, as is partly the case for N = 
theories |J. 

In a preceding paper we addressed these questions in the simplest cas e k = 1. We showed 
that the best point of view was to consider twisted N — 2 Q, [t2| , and that realizations of such 
theories were nicely provided by the geometrical problems of polymers and percolation. N = 2 
supersymmetry was manifest in the structure of the correlators and the geometrical operators 
algebra. The presence of supersymmetry was traced back to the Parisi Sourlas |l3fl supersymmetry 
of the lagrangians (involving bosons and fermions) that are used in the description of geometrical 
models (n — > limits). The knowledge of supersymmetry allowed us also to predict some new 
exponents, like the fractal dimension of the backbone of percolation in two dimensions. 

The purpose of the present paper is to extend these observations to the case of arbitrary k. We 
find in particular lattice analogs of the chiral ring in integrable geometrical models, and families 
of multicritical polymer points connected by the N = 2 flow. 

The organization is as follows. In the second section we work out the reexpression of N = 2 
partition functions in terms of generalized Coulomb gas that was started in |5|. In section three 
we study the T)~ vertex models using a new geometrical formulation that generalizes the k = 1 
case. This formulation involves bound states of strands carried by the edges of the square lattice. 
It exhibits very clearly the symmetry properties of the q — exp^in/k + 2) points, at which the 
continuum limit is the Ak+i modular invariant for the N — 2 series. For instance the free energy 
can be obtained in a simple way at that point, without using any Bethe ansatz computation. The 
continuum limit is indeed derived. We also identify a set of lattice operators that reproduce the 
chiral ring in an appropriate continuum limit. In the fourth section, we consider the duals (in 
the sense of Coulomb gas in two dimensions) of these models, and argue that they correspond 
to multicritical polymers. Generalizing the analysis of fl0|| , the various sectors and the role of 
fermions and bosons, are geometrically interpreted. We discuss the coincidence of the exponents 
formulas with the so called Flory formulas u% . 



3 The N = and N = 1 minimal series also possess in principle Landau Ginzburg descriptions, however these 
cannot really be used to obtain precise information on the system, like critical exponents. 
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2 Generalized Coulomb Gas Expressions for the N = 2 Par- 
tition Functions 

The generalized lattice Coulomb gas was introduced in |)j to describe the continuum limit of SU{2) 
minimal coset models and of Tk vertex models based on spin j — k/2 representations of U q sl{2). 
It involves as expected |fl5f a free bosonic field and Zk parafermions, which are coupled through 
boundary conditions. Introduce the partition function of a free bosonic field with coupling g and 
windings 2mir, 2m'ir along the generators wi,u>2 of the torus, with r = UJ2/0J1 



Zmm<{g) = \l ~r- — =exp (--^— |m - to't| 2 ) (1) 
V i itit r\r\ V 1 rrir / 

where rj is Dedekind's function. Introduce also for a simply laced Lie algebra G with Coxeter 
number k+1 the partition functions Z^(r, s) of the Zk model of |l7j with boundary conditions 
twisted by e 2l7rr / fc , e 2t7TS / k ( W e will not need their exact expression in the following). The generalized 
Coulomb partition function was then defined in || 

Zg(g)= ]T Z%{r,s) Z mm ,{g) (2) 

r,s=o,...,fc— 1 m = r mod k 

m' = s mod k 

with central charge 

«-ST5 (3) 

In the following we shall concentrate on Lie algebras of A type, and henceforth suppress the 
superscript G. The coulomb partition functions satisfy the duality relation 

Z gc {g) = Z gc (^ (4) 

It was then shown in that the N = 2 partition functions for the minimal Ak+i series with 
central charge ^ obtained in |l8| could be rewritten as follows^ 

Z'N=2.k — 2 (Zap + Zpa + Zaa) 

fl \ fc + 2 

2~2 Z k(r,s) +S mAm , mod 2 J Z mm ,{-^-) (5) 

r,s=0,...,fc-l m _ r m0( J fa 

m! = s mod k 



The doubly periodic contribution is also known Jljto be 

Zp P = Tr{-) F = k + 1 (6) 

4 The labels A (antiperiodic) and P (periodic) refer to the fermion boundary conditions on the torus 
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Notice that this index does not vanish. 

We can now simplify the expression (|^) using the following arguments. It was also shown in 
|| how the generalized Coulomb gas can be used to reexpress the partition functions of minimal 
models 

SU(2) k ®SU{2)i- 2 (] 
SU(2) k+l _ 2 ( ' } 

with central charges 

3k 6k 



k + 2 1(1 + k) 

This reexpression involves the generalized Coulomb gas with coupling 

a = \-TTk (9) 

plus some interactions between the winding numbers. Consider now the particular case I = 2 that 
is associated with lattice models of type Ak+i- For such models c = and one expects the partition 
function to be a constant, equal to k + 1 (the number of possible heights per lattice site). On the 
other hand, following || one gets 

z A k+1 = 12 Zfe ( r ' s ) 12 Zmm ' ( fc(fc 2 + 2) ) 12 cos { 27r kT2 m A m ') (10) 

r,s=0,...,fe-i m = r mod k n=1 

m! — s mod k 

The sum over exponents produces a (k + 2)S mAm , moc [ fe+2 — 1- Analysis of each of the two terms 
and use of the duality symmetry gives therefore the formula 

Z Ak+1 =k + l= Z ^ s ) E (2^ Aro 'niod 2 -l)^wf^ N ) (11) 



o,...,fe l m — r mod k 

m! = s mod k 



Notice that this can be reformulated as 

\ ' ( _ynArn' Jr __ j ( [ 

2k 



k + l= 12 Z k (r,s) £ (-r Am 'Zw(^r) (12) 



r,s— 0,...,k l m = r mod k 

m' = s mod k 



Z' N=2 , k = Z gc [ — )+ — = Z gc [ Tn - T ^ ) + — (13) 



This identity generalizes the Euler's identities. One deduces from (||) and ( |ll| ) 

'fc + 2\ fc + 1 / 2 \ . k + 1 

~2k~) + ^2~ ~ 9C \k{k + 2) 

The full N = 2 partition functions with projection on odd fcrmion number writes therefore 

'k + 2 
~2~k~ J ~ ~ 9C \k(k + 2) 



Z N=2 ,k = Z gc [ — )=Z gc [ T71 - T ^ ) (14) 
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This identity was originally derived in p9| and independently in pC[ | . It provides a first clue con- 
cerning the lattice models whose continuum limit is N = 2 supcrsymmctric. Indeed the generalized 
Coulomb partition functions are the continuum partition functions for the Tk vertex models based 
on the quantum algebra U q sl(2) and spin j — k/2. However this identification leaves totally ob- 
scure the origin of the N — 2 supersymmetry. At first sight, nothing distinguishes the T& vertex 
model at the value q = exp(iir/k + 2) from any other point on its critical line jH). The next section 
is devoted to a geome trical formulation that will clarify this issue. 



3 Geometrical Reformulation of the Tk Vertex Models 
3.1 Generalities 

The weights of the Tj, vertex model are conveniently encoded in the R matrix. This operator acts 
in a space (C fc+1 )' 8)2 , and its elements are the Boltzmann weights of the corresponding configura- 
tions (figure 1). An important technical point is that several choices of weights are possible that 
are obtained by gauge transformations. Such transformations keep the Yang Baxter equation 
satisfied. Moreover they do not affect the torus partition function of the system. There is a special 
transformation that ensures commutation of the R matrix with the generators of U q sl(2), and in 
the follow ing we always consider this R matrix as giving the weights of the model. When the R 
matrix is in the commutant of U q sl(2), it satisfies some nice a lgebraic properties that define for 
instance in the j = 1/2 case the Temperley Lieb algebra. An advantage of focusing on these 
algebraic properties is that other representations can then be used, which make the physics of the 
problem more transparent. 

Consider for instance the j = 1/2 case. The corresponding 6 vertex model has a R matrix 
(called r in the following) that writes 

r = I + f(u)e 



where 



u is the spectral parameter, and 



exp(«7),7 S [o,7r] 



/(«) 



sin(7 — u) 

The operator e writes in the usual spin j = 1/2 basis 



(15) 



(16) 



(17) 



/ 





°\ 


q- 1 


-1 





-1 


9 





{o 





o J 



(18) 



and satisfies the equations 



(q + q- 1 ) 



0,\i-j\ >1 



(19) 
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that define the Temperley Lieb algebra. In these equations, e% is meant for the operator that acts 
as ( p"g| ) in the i th and i + 1 th copies of C 2 , and identity otherwise. 

It is however convenient to think of another representation provided independently by knot 
theory ideas pi] , and by polymer and percolation problems. We notice that the Temperley Lieb 
algebra can be represented graphically as in figure 2. The operators act now in a space of lines, the 
identity just propagates the lines straight, while insertion of the Temperley Lieb matrix e attaches 
two neighbouring incoming and outgoing lines. Every time a closed loop is formed it is assigned 
the value f\ [2] q = q + q~ 1 . We can therefore think of the 6 vertex model as a loop model. There 
are now only two local configurations, represented in figure 3, which we call (this notation will be 
explained later), I, x, with respective weights 1, /(it). A very important remark is that on a torus, 
the partition function of the loop model equals that of the original 6 vertex model if and only if 
the non contractible loops get a weight 2 instead of [2] q . 

The price to pay for having now only two local vertices is that the computation of the weight 
of a configuration requires in general a non local information, namely, when do pieces of lines 
attach together to form a loop, that should then be assigned a weight [2] q . We recall that the loop 
evaluation used here is the geometrical equivalent of the Markov trace p2[ . There is however a 
magic point where this non local information is not needed, it is obviously when [2] q = 1. This 
corresponds to q = exp(iir/3), g = 1 — 1/3 = 2/3, ie it is the appropriate value for the N = 2 point 
on the 6 vertex (Gaussian) line. Notice that at that point, the free energy can be evaluated in a 
trivial way: every vertex has two possible states, so in the large volume limit where we can neglect 
the special weight of non contractible loops we have 

/ = log 2 (20) 

This point was actually studied in jl0| from the point of view of twisted N — 2, and connected 
with the percolation problem. However, it does not seem that one can consider the other values of 
k in the series as representing some multicritical percolation points. The situation will be different 
for polymers, see below. 

3.2 Symmetrizers and Fusion 

The Lfe vertex model can be obtained from the 6 vertex model by projecting k spins j = 1/2 onto a 
spin k representation of U q sl(2) p3|, p4[. T he necessary object for doing that is the q symmetrizer, 
which has been widely studied p4 |26|, |27[ . It can be expressed using Temperley Lieb operators as 
follows. First recall that 

Po 



1 



q + q- 1 



Pi = S 2 = l-P (21) 

and introduce 

s = q- l Px- qPo (22) 



5 In the following, [x] q denotes as usual the q-analog, equal to (q x — q x )/(q — q x ) 
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Then one has 

[k] q 



-fc(fc-l)/2 



where the sum runs over permutations of k objects. Decomposing a as a product of transpositions 
Ti t i-\-\ gives the set of indices I a , of minimal total number \I a \. The symmetrizers satisfy some simple 
induction relations that make their determination easy. Now notice that S 2 can be represented 
graphically ]2£| as in figure 4. The property PqP\ = translates into the fact that strands always 
have to pass through £2, and cannot make a U turn. This property generalizes to any k. 

Now the fusion procedure for getting the Boltzmann weights of say the T 2 = 44 vertex model 
can be implemented in a variety of ways. We follow here the analysis of |^5|, ^7) . One finds in these 
references the proof that 

Rj=i = 5 2 <S 2 r 2 (u - ~f)r 3 (u)ri(u)r 2 (u + j)S 2 S 2 (24) 

This a natural generalization of the well known formula used in conformal field theories or knot 
theory, that holds in the u — > ioo limit. It corresponds to the picture of figure 5. Formula ( pi[ ) 
generalizes easily to give the R matrix of higher k vertex models as well. One finds 

Rj=k/2 = S k S k {r k [u - (k - 1)7] . ..r 2 k-i[u}} x {r k -i[u - (k - 2)7] . . .r 2k _ 2 [u + 7]} 

x ...x{n[u]...rk[u+(k-l)i\}SkS k (25) 
We remind the reader that R has also an expression in terms of the quantum projectors p9| 

R ( u ) = p ^ + r^^- 1 + • • • + i-^i-^- 1 ) • ■ ■ t^v Pj=0 (26) 

where y = exp(-iu). 



3.3 Geometrical Representation of the Tk Vertex Models 

Let us start by discussing again the j — 1 T 2 = AA vertex model. We consider the strand repre- 
sentation of the Temperley Lieb algebra introduced above, and follow it through fusion. First it 
is clear that the R matrix acts now on bound states ie on states that have two lines incoming 
each leg of the four vertex. The symmetrizer acts on these lines, which is indicated by a transverse 
bar. Now we can compute the product (|24|) by picking one of the two terms in each of the r 
factors, and drawing accordingly either lines going straight or making some U turn. There are 
many allowed such configurations, that are represented in figure 6rl However after acting with the 
symmetrizers, only 3 configurations remain, which we call 1 figure 7. Their respective 

weights are W\, W x , W x 2 , given by 

W 1 = 1 

W 2 = /( U -7)/(w + 7)/ 2 W 
W 3 = fju-^ + f ju + 7) + 2f(u - 7 )/(u + 7) [f(u) + cosy] (27) 

6 This computation appears originally in some unpublished notes of Lou Kauffman devoted to "Multi strands 
calculus" and cabling in knot theory. The present formulation with insertion of symmetrizers leads to graphical 
rules for the SU(2) generalizations of the Jones polynomial (L. Kauffman, H.Saleur, work in progress) 
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The isotropic point is, as for any k, obtained for u = 7/2, and one has then 



Wi = 1 

W 2 = 

l + [2], 

W 3 = 1 (28) 

so that the W are invariant by a 90° rotation. The three allowed configurations can be written in 
the Temperley Lieb algebra as 

1 <-> S 2 S 2 
x «-► S2S 2 e2S 2 S2 
x 2 «-> S 2 S2e2eie 3 e2S2S2 (29) 

As for the 6 vertex model, although the interactions are much simpler in the geometric picture 
(44 vertices replaced by 3), the computation of the total Boltzmann weight of a configuration 
involves complicated evaluation of loops weights. In our case, a single contractible loop has weight 
[2] q , while a symmetrized contractible double loop (see figure 8) has weight [3] q . A single non 
contractible loop still has weight 2. The value of q at which one can forget the loop weights, and 
the geometrical model becomes localQ is given by 

[3],= 1, q = exp(iir/A) (30) 

There are various ways of proving this assertion. One is to argue that since the q dimension of the 
spin one representation, which is the fundamental representation to consider in this case, is equal 
to one, the full Hilbert space onto which the transfer matrix is acting has also a q dimension one, 
irrespective of the number of spin one representations that are considered. Therefore all quantum 
symmetric observables act trivially onto this space, and there cannot be any non local contribution. 
This can be made more precise, but is not very geometrical. Another proof is elementary and uses 
the very geometrical definition of the model. Consider a typical "fat" graph appearing in the 
computation of the partition function of a large system (see figure 9). Isolate in this graph an 
external "ear", and expand the symmetrizers that sit on the two branches. One gets the pictorial 
equation represented in figure 10. The first term has prefactor 



2 [3], - 1 



[2], [2] 



v 



that vanishes precisely at the point q — exp(i7r/4) = exp(z7r/ (k + 2)). At this point we can therefore 
simply cut each ear to evaluate the value of the graph, counting a factor l/[2]^ for each suppressed 
ear. By repeating inductively this procedure we end with a graph that has no ear and therefore 
looks as the one of figure 11. We can decimate the loops using the pictorial equation of figure 12 

7 The model is truly local for an infinite system only. In a finite geometry, care has to be taken with non 
contractible loops. They however , for the present regime, are negligible for the evaluation of thermodynamic 
quantitites 
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to get a factor [2] q — l/[2] g per loop, , which equals l/[2]? in that present case. The last loop gets 
a factor [3] g . Therefore the total weight of the graph turns out to be 



I \ area-i 



™graph= (j^J ( 31 ) 
Now notice that 

area = number of vertices of type x + 1 (32) 

Therefore we can evaluate the partition function by forgetting about the weights of contractible 
loops and multiplying the weights of each vertex by a factor 

Ai = 1, A 2 = l/[2]q, A 3 = 1 (33) 

It is also interesting to consider the correspondence between the geometrical vertices and the 
expansion of the R matrix in terms of projectors. One finds 

1 «-> Pj=0 + Pj=l + Pj=2 

i 2 i«-i4K i+ ( [21 '-i4K° 

x 2 « (Pi; -1)^.0 (34) 

The extension of these results to higher values of k is easy. In general, the legs of the vertices 
carry bound states made of k symmetrized lines. There remains only k + 1 vertices, that can be 
labelled l,x, . . . ,x k see figure 13. These vertices correspond to interactions that are once again 
more conveniently written graphically. We draw as in figure 14 a square with each face carrying 
an e matrix with the appropriate label, and multiply the e's from top to bottom, and inside a 
given line, from left to right. The precise form of the weights is not needed here (it involves use 
of quantum 6j symbols [^0)). They have the property of being symmetric under 90° rotation for 
u = 7/2. At his point one has therefore the symmetry 

W x i = W x *-i (35) 

The symmetrized contractible k loops acquire a weight [k + l] q that becomes equal to one for 

q — exp(iir/k + 2) (36) 

At this point one can show that the geometrical model becomes local, with weights renormalized 



by the factors) 30 



A,. = ^-y- (37) 

Notice that at the special point ([36j) one has 

A x » = \ x k~i (38) 

Notice that, as in the k = 1 case, the free energy can be simply evaluated at the N = 2 point of 
the Tk vertex model. Since each vertex can be in one of the k + 1 allowed states independently 
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of its neighbours, one gets in the large volume limit where the effect of non contractible loops is 
negligible 

f = log[J2K*W x A (39) 



3.4 Torus partition functions 

The loop model introduced above being equivalent to the Tk vertex model without charges at infin- 
ity, it is easy to derive its continuum limit using the generalized Coulomg gas mapping. Following 
the results of || one has 

9=1 k-^=W^2) (40) 

The coupling between the parafermions and free bosonic field quantum numbers is built in in the 
model. The only remaining point concerns the allowed range of to, to'. Because we want the x 1 
to be true operators acting in a transfer matrix formalism, we need the square lattice at hand to 
be described diagonally as in figure 15. There is therefore an even number of edges crossed when 
describing one of the generators ct>i,£d2, so that following ]f|, to, to' are integers. Formula ( |l4| ) 
therefore follows. 



3.5 Lattice Candidates For Chiral Primaries 

We now present some arguments that the x % operators, that become local operators at the magic 
point q = exp(iir/k + 2) for the j — fc/2, k strands case, are lattice candidates for the chiral 
primaries |pl| . Of course as lattice operators, the x 1 have a left and a right moving part, that 
behave identically. It is more proper to identify their lattice algebra with the direct product of 
left and right chiral rings of the continuum limit. For simplicity we think of the left moving sector 
only in the following. 

It is straightforward to perform the multiplication x l .x 3 . For this one simply connects all the 
outgoing lines of a x J vertex to the incoming lines of an x l one, and expands the intermediate 
symmetrizers (figure 16). The result takes the form 

min(i+j,k) 

x\x 3 = °i xl ( 41 ) 

1=0 

The precise form of the q coefficients is not needed in the following. The important point is that 
the right member truncates and does not contain terms with exponents greater than i + j. The 
reason for this is not totally obvious. To form a term x l requires the introduction of I matrices 
efc (figure 17). On the other hand, x % and x J contain respectively i and j such terms, while the 
symmetrizers do not contain any. Therefore there are not enough terms to form an x l when 
I > i + j in the product. The fact that i + j cannot exceed k is simply because we start with k 
strands on each leg, so we cannot connect more than 2k of them. 

Now we can try to transform such a product into a meaningfull expression in the continuum 
limit. For convenience we suppose that we work on a cylinder with some strands states propagating 



9 



in the time direction (figure 18). After the conformal mapping from the plane to the torus, fields 
become operators acting on the strands. Due to the conformal factors, these operators have matrix 
elements that scale like L~ h ~ h , where L is the radius of the cylinder. On the other hand the discrete 
operators x l introduced above have matrix elements that scale like 1 . Therefore we define rescaled 
operators 

X l L = L- 2h ^x i (42) 

so the product ( pil"| ) reads now 

i+j 

X l L .X{ =Y J CiL^ hi +i- l - hi ~ h ^X l L (43) 

1=0 

This equation has a (non trivial) limit when L becomes large if and only if 

h i+j = hi + hj (44) 

In that case fli^ ) becomes in the large L limit 

X\X j = X l+J for i + j<k, otherwise (45) 

This algebra is precisely the chiral ring 0, || , while the condition ( f44| ) occurs as well in the study of 
continuum N — 2 theories when one requires the short distance expansion to reproduce the chiral 
ring. One finds then 

h x i = i/2(k + 2) (46) 

Another point concerns the U(l) charge of the x % operators. Following |l(| let us define it as 
proportional to the maximum number of non contractible loops (in the time direction) the insertion 
of x l can create. One finds 

Q xi cx i (47) 
again a result known to hold for N = 2 theories where 

Q x <=i/(k + 2) (48) 

It is as well possible to consider the x l as antichiral primaries. They satisfy the same algebra, 
but we define their charge as proportional to minus the maximum number of non contractible 
loops (in the time direction) their insertion can make contractible (figure 19). Notice also that 
the 90° rotation formally correspond to the so called Poincare duality [|[ Finally it is clear 
that I behaves indeed exactly as the identity in the transfer matrix formalism, and therefore has 
vanishing conformal weight. We discuss the conformal weights of the other x % fields in the next 
section. 

3.6 Discussion 

First we discuss the case j = 1/2. We therefore have contractible contours with weight [2) q = 1 
for q = expi7r/3, non contractible contours with weight two, and two vertices of interaction, I,x 
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with the same weight one. We consider the model on a cylinder. It is easy to isolate the doubly 
periodic sector by considering for a while configurations where the non contractible loops also have 
a weight one. In that case the partition function is obviously equal to a constant, in agreement 
with Zpp — = 2. This last value is merely an overall normalization problem when taking 
the continuum limit in the Coulomb gas mapping. It can be found by noticing that the loop 
configurations so isolated just compute the partition function of the A 2 model, which corresponds 
in turn to taking yet another representation of the Temperley Lieb algebra of solid on solid type, 
with elements for the face configuration of figure 20 



[h\q 

We can therefore picture the Ramond ground state as a frustration line running accross the cylinder 
in time direction, such that every loop crossing it gets a weight 1 instead of 2. In the Coulomb gas 
mapping, this corresponds to having a pair of electric charges at the extermities of the cylinder, 
with value eo = ±1/3. This has dimension h = e^/Ag = 1/24 as expected. Let us now consider the 
full theory, with non contractible loops having weight 2^ . Because insertion of I acts identically 
on the propagating strands, its dimension is certainly 

hi = (50) 

Consider now the operator x. It is reasonable to assume that the connected correlation function 
with its antichiral conjugate is obtained by creating a pair of non contractible loops where x 
is inserted, and destroying it when the conjugate is inserted: in other words by selecting the 
configurations where a loop connects the two insertions as in figure 21. It is well known how to 
evaluate this dimension in the Coulomb gas mapping fl32| . The corresponding operator is magnetic 
with charge m = 1, and therefore 

= 9 {m=lf = 1 

4 6 v ' 

as expected. This dimension is also what would be observed if we were to compute the correlation 
function in the antiperiodic (NS) sector. 

This picture essentially generalizes to the higher k case. The doubly periodic sector is repro- 
duced by using the restricted solid on solid model representation of the Temperley Lieb algebra 
for The Witten index is k + 1 corresponding to the number of nodes of the diagram. The 

Ramond ground state is obtained by giving to the non contractible loops accross the cylinder a 
weight [2] q . Finally if we admit the likely hypothesis that connected correlation functions of x l 
operators and their conjugates are described by i loops joigning the two insertions we find, using 
the generalized Coulomb gas mapping, the weights h x i = i/2(k + 2). 

4 D type invariants 

It is well known how the D 2+k / 2 theories can be obtained from the A^+i theories by a Z 2 orbifold. 
In our case it is easy to form symmetric combinations of the diagrams x % and x k ~ l . However we 

8 Unfortunately we do not know what is the microscopic definition of all the sectors in that case. In that respect, 
the polymer situation described in the next section is more favorable 
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do not know what the equivalent of the y field [g, g| should be. Hints in the identification of D 
lattice models can be obtained from consideration of the partition functions. Using once again the 
results of one finds first that the expression (Q) generalizes to D or E invariants, but with 
the appropriate parafermions partition functions (0). Then since D parafermions are themshelves 
orbifold of the A ones one has 

Z° 2+k/2 (r, S ) = ±[Z k (r, S ) + (~yZ k (r,s + k/2) + (~yZ k (r + k/2, S ) 

+ (-) r+s Z k (r + k/2,s + k/2)] (52) 

The peculiar combination of sectors of the T k vertex model reproducing the D 2 + k /2 invariant can 
obviously be read from the above equation. But it still lacks a natural interpretation. 

5 K th Critical Polymer Models and Twisted N = 2 with 
Level k 

There is another, and physically more interesting, way of interpreting the N = 2 partition functions. 
Recall that the T k vertex models are expected to have in fact two phases. In the first one, the 
model is integrable with Boltzmann weights given by the usual formulas p3| , p4| , and coupling 
constant as above 

9 = \-1, 9 = e ' i7 : 7£ [0,1/fe] (53) 

The second phase is not very well known except for k = 1 . It is supposed to be obtained by adding 
some vacancies. There are however some strong arguments that the corresponding renormalized 
coupling constant is 

.9=^+7 (54) 

In the case k = 1 the first regime corresponds to the critical Potts model with Q = (q + q^ 1 ) 2 , 
the second regime to the critical 0(n) model with n = q + q^ 1 . For higher k it is likely that 
the second regime de scribes a sort of multicritical 0(n) model with in particular the formation 
of boun d states, ie up to k lines colliding on a same edge, allowed. Now choose 7 = 1/2. This 
value corresponds to n = ic to k th critical polymers. The associated coupling constant is 

k + 2 

g = "IF (55) 

dual to the coupling constant for the vertex model that we discussed in the preceding section. In 
the case of polymers, we can however generalize the discussion of [To| for the dense and critical 
case, and provide a natural geometrical interpretation of the various quantities so far encountered. 
The natural point of view here is to consider twisted N — 2 theories JT^ . Recall that twisting 
the theory is obtained by adding a term |9 J to the stress energy tensor. The new central charge 
vanishes. The partition function of the twisted model is however the same as the one of the 
untwisted model, so twisting merely amounts to a change of point of view. The former ground 
state of the Ramond sector is now considered as the true SL2 invariant ground state, while the 
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former true ground state is now considered as a state with negative dimension. Notice also that, 
because of twisting, the fermio ns acquire integer dimensions and therefore have identical boundary 
conditions for the cylinder and the plane. 



5.1 Polymers Partition Functions 

First we notice that for polymers also the bulk free energy is easy to evaluate. Polymers occupying 
a vanishing fraction of the available space one has / = 0, as is usual for n — > limit models. 

As in jL0| we think of the n — ► limit as obtained by associating to each polymer loop a 
bosonic or a fcrmionic variable. In this way all contractible loops on a torus disappear, while 
non contractible loops can get a weight zero or two depending on their winding and the fermions 
boundary conditions (figure 22). One has, if the labels A or P refer to these discrete fe rmions 
boundary conditions 

Definition 

^AP(resp.PA,resp.AA) — Sum over configurations with an even number of non contractible 
k th critical polymers of total length C, that cross ui^ (resp.wi,resp.a;i + 0J2) an odd number of 
times, with weight 2 mxm} ° er of Polymers^-£ x t j ie Boltzmann weights induced by the k th critical 
interactions 

where ijl is the appropriate inverse connectivity constant. One can also define the manifestly 
modular invariant quantity 

Z e = I Sum over configurations with an even number of non contractible k th critical 

polymers of total length C with weight 2 number of P°lymers 
fi~ c x the Boltzmann weights induced by the fc tft, critical interactions} (56) 
One checks easily from these definitions thatQ 

Z e = - (Z A p + Zpa + Z A a — Z PP ) (57) 

Notice that projection on odd fermion number occurs naturally on geometrical grounds. 

On the other hand the analysis with the generalized lattice Coulomb gas is such that topological 
properties (magnetic defects, electric charges, winding numbers...) do not depend on the peculiar 
k th critical regime. We can therefore use the same expressions as the ones in to express in the 
continuum limit 

Zpp = Z k (r,s) ^2 Z mm i 

r,s=0,...,fc-l m — r moc |/j 



1)1 



= s mod k 



z AP = £ z k (r,s) Yl (-r(-r Am 'z n 

m = r modfc 
ml = s mod k 



r,s=0,...,fc-l 



3 Recall that the two winding numbers of a non contractible loop are coprimes 
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z PA = y, Zk ^ s ) E (-r'(-r Am 'z n 

m = r modfc 



r,s=0,...,fc-l 



Zaa = Z ^ s ) E H^" 1 ' (-) mW Z mm , (58) 

r,s=0,...,k-l m = r mod/c 

w! = s mod k 

One checks that each of these expressions coincides respectively with the known partition functions 
of the N = 2 theories. In that case we therefore have found a geometrical interpretation for all the 
sectors of the theory. Let us emphasize again that due to twisting, the fermions acquire integer 
dimensions and therefore have the same boundary conditions in the plane and on the torus^. 
Summing these four contributions one gets 

( k + 2\ 

Z e = Z gc ( ^ J , for k th critical polymers (59) 

As in the k — 1 case, this partition function is not consistent physically because it does not contain 
the identity (for the twisted theory), as expected since we projected on odd fermion number. The 
physical polymer partition function is rather obtained by projection on even fermion number 

Z phy = Z e + k+l^Z e + k + l (60) 

The ground states in the twisted theory describe simply configurations with no polymers. That 
they appear with some multiplicity occurs probably because of the additional degrees of freedom 
introduced to make the polymers multicritical (vacancies and their generalizations J33|). As in 
the k — 1 case we notice that R and NS correctly describe observables with an even number 
of polymers only. The full polymer theory would need i n addition the consideration of a sector 
with Z4 twists |l(J , which is easily studied via the spectral flow. Also we recall that the BRS 
cohomology used to obtain topological theories out of twisted N — 2 jjjj is not the right procedure 
to extract physical states from the polymer point of view. On the contrary, Q bus turns out to be 
the operator that creates polymers out of the vacuum. 



5.2 Some Polymer Exponents and the Flory Formula 

The most important quantity is the exponent v that controls the mean size of polymers. Consid- 
ering a polymer loop of C monomers one has asymptotically 

<i£>oc£ 2 " (61) 

where Rq is the radius of gyration. If we introduce the 2 legs polymer operator, of conformal 
weights h = h, whose correlations are defined by summing over configurations of a polymer loop 
attached at two different points, one has, by standard scaling arguments 

- =2-2/i (62) 
v 

10 In the preceding paper IICJ] we referred to R and NS as the boundary conditions in the plane 
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This operator was identified in pToj j to be the first non trivial operator in the periodic sector (X in 
the Landau Ginzburg picture), with weight after twisting given by 

h =kh ^ 



Using formula ( |62[) one finds 

k + 2 



- a( * + i) (64 » 

These exponents have already appeared in the literature and they are known as the Flory exponents 
Jlij for k th multicritical polymers (in our conventions, ordinary polymers are k = 1). The Flory 
formula is obtained by making crude, and notably wrong [ p4| , assumptions for the free energy 
of a polymer. One usually writes F as the sum of an elastic and energetic contributions, and 
minimizes with respect to say C at Rq fixed to obtain v. The elastic contribution is argued to be 
F e i (x Rq/£, and the energy contribution for a k multicritical point F en cx Rq(C/Rq) +1 where 
d is the space dimension |35| . The argument for this last expression mimics field theory where for 
ordinary polymers, the (/> 4 term dominates, corresponding to two pieces of polymers that come in 
contact (figure 23), for tricritical polymers the cf> e term dominates corresponding to three pieces of 
polymers coming in contact (^6)... Writing therefore 



Rn ~J ( £■ 



fe+1 



one finds after minimization 



f «^ +r Km) (65) 



k + 2 



~ 2 + dk (66) 
One could not really take this formula seriously based only on the Flory " derivation" . However 
in the case k = 1 it is in fact exact for d = 1,2,4 and very close to the numerical estimates 
for d — 3. v = 1/2 for d = 4 which one can therefore correctly identify as the upper critical 
dimension. Formula ( |66[ ) should not be applied above the critical dimension, where the exponent v 
is expected to remain equal to 1/2. Some explanations of approximate validity of Flory's formula 
have been proposed j37j]. What is fascinating however is that it can be exact in most cases. So 
far higher k's had not been really considered. We notice however that if we put d = 2 in (^) we 
recover the exponents (^34|). This shows they are also exact for the peculiar kind of multicritical 
polymers at hand, and strongly suggests that the Flory formula can be exact because of non 
renormalization theorems due to the hidden N = 2 supersymmetry in polymers. As a final 
argument in that direction, notice that v in ( |66|) becomes equal to 1/2 for 

d = 2 + \ (67) 
k 

which is precisely the upper critical dimension for k th critical N — 2 Landau Ginzburg theories. 
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5.3 Renormalization Group Flow 

It is worth considering the problem from the point of view of renormalization group flow |], [IJ. 
We remark that because twisting is only a change of point of view on the system, the flow results 
established for the untwisted case should hold also for polymers. Indeed perturbing for instance by 
X the k theory means shifting the geometrical weight of the polymer to (n + Sfi)~ c , independently 
of whether we consider the ground state as the state without polymers, or the state with large non 
contractible loops that have weight two. 

For the N = 2 series, the "minimal" supersymmetry preserving flow is a flow from the theory 
with fc to the next theory with k — 2 (X k ~ 1 is redundant). It is generated by the superpartner of 
the top component of the chiral primary field X k , with expression in the generalized Coulomb gas, 
where ipi is the parafermionic field of lowest dimension, 

$ pert = Tpltfiexp {j^ z ' J ^j ( 68 ) 

and dimension h = h = After twisting, since G + acquires weight one, this field has the same 
dimension as X k in the twisted theory, which can also be computed by adding the proper charge 
at infinity for the field 4> m the generalized Coulomb gas. One finds 

This perturbation is still integrable in the twisted theory as discussed in [[l2| . The above dimension 
provides the value of the exponent v u for the multicritical polymers at hand [ |43[ 

vu = ^ (70) 



and therefore the crossover exponent 



tp = — = -T-— (71) 
v u k + l 



We recall that the renormalization group flow in the N = 2 series is more subtle than in the N = 
or N — 1 series. The analysis of the superpotential X k+2 + \X k can easily be done, at least 
naively, since there is only wave function renormalization plfl . However two points with different 
values of k are infinitely far away in the sense of |l2| , so the flow cannot be studied in the conformal 
perturbation framework^. This is because the index takes value k+l for the k th critical polymers, 
and therefore has to jump by two units between the k and the fc — 2 theory. As long as < k < oo 
one does not expect spontaneous breaking of supersymmetry. 

The perturbation with the operator X is also integrable Jf4| , and physically interesting. Recall 
that X was identified with the two polymer legs operator, and therefore adding an X perturbation 
is like changing the geometrical weight of polymers. Two behaviours can occur depending on the 

11 From the polymer point of view, and owing to experimental and numerical knowledge of multicritical polymer 
systems, there is little doubt that such flow occurs 
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sign of the coupling. For the sign corresponding to <5/j > 0, the theory flows to a trivial fixed point. 
This describes so small polymers that they disappear at large scale. For the sign corresponding to 
<5/j < 0, we expect physically the polymers to become so large that they occupy a finite fraction 
of the available volume, with fractal dimension exactly equal to two, that is an exponent v = 1/2. 
This means the system should flow to the dense polymers phase, that is described by an 77, £ system 
fl(i|| , and that supersymmetry is spontaneously broken. The decoupling of the ground state has 
a very transparent physical interpretation. Indeed for the state without polymers Z oc 1, while as 
soon as one dense polymer is allowed Z oc e / xarea for f some non vanishing free energy. Notice 
that in the dense phase, the free energy is now a non trivial number, which is not known exactly 
except in special cases (hamiltonians walks on the Manhattan lattice for instance). 

For k — 1 (ordinary self avoiding polymers), the preturbations X k and X coincide. The 
"bottom" of the series of multicritical points is therefore dense polymers, where the supersymmetry 
is spontaneously broken. 



5.4 k — > 00 limit 

Another insight can be obtained by considering the k — > 00 limit. In that limit the exponent v tends 
to 1/2, the exponent of dense polymers^} This corresponds to polymers that are collapsed onto 
themshelves, due to strong attractive interactions, or compression from the exterior. We therefore 
find again dense polymers, which are indeed expected to be the end point of the multicritical 
polymers series when described in the direct renormalization or 0(n),n — > multicritical field 
theory. We can here make the identification of the k — > 00 limit of our theories with dense 
polymers very precise. Indeed it was shown in pX[ that dense polymers are conveniently described 
by an 77, £ system with central charge c = —2. This corresponds to breaking the symmetry between 
bosons and fermions in say the free field representation of N = 2 twisted theories. By summing 
over various sectors, it was shown that the entire partition function of dense polymers with an 
even number of non contractible loops on the torus is a gaussian partition function 

Z dense = Z cm (72) 

where by Z c we mean an expression similar to Z gc but with the partition functions of the Z\ 
parafermions taken to be equal to one. On the other hand it is possible to work out the k — > 00 



limit of the partition functions (14). For this purpose notice first that the coupling constant goes 
to a finite value g = 1/2 in that limit. For k large enough, the 771,777/ excitations are rapidly 
damped out in the second summation in the expression (|2[) of the generalized Coulomb partition 
function, so we can truncate this summation to m = r, m = s. Moreover as k becomes large, the 
frustrations of the Zk model become negligible for the values of r, s such that Z rs is not itself too 
small. One can therefore write 

Z gc [^j^j « Z k (0,0)Z c [l/2], as k 00 (73) 

Now such a result can also be established for the level k Wess Zumino partition functions which 
can be rewritten as a generalized Coulomb partition function for a coupling g = 1/k. On the other 

12 u = 1/2 is of course also the exponent of brownian walks, but this is not apparently the same problem, although 
some other features are common with the dense case. 
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hand, using reexpression of these same WZW partition functions in terms of 3 bosonic partition 
functions as worked out in Q leads to 

Zwzw oc ( / -}— _ ) (74) 
V v IraT-qri J 



Combining these results we get finally 

Z 9c J oc ( rr— _ ) Z C [V2] as k -> oc (75) 



fc + 2 N ' ' ■ 2 



2k J \ %J IniTr/f] 

Therefore in the large k limit we find indeed the dense polymer partition function, up to a free 
bosonic field factor that helps maintaining central charge at a value zero. As in the preceding study 
of the X perturbation, one expects in the k — > oo that spontaneous braking of supersymmetry 
will occur, leaving only the ij, £ system with partition function Z c [l/2]. Notice that the untwisted 
model corresponding to the k — > oo limit is the c = 3 N — 2 theory, which is well known to 
have vanishing index. Therefore dense polymers occur both at the "top" and the "bottom" of the 
multicritical series. That dense polymers appear both in the k — > oo and in the low temperature 
phase is due to the existence of two independent mechanisms for obtaining them. When k — > oo, 
more and more attractive interactions are added, so the polymer is forced to collapse onto itself for 
internal reasons. On the other hand, the breaking of symmetry of the 0(n),n — ► model actually 
reduces the available volume, so the polymer collapses for topological reasons ffof . 

5.5 Ordinary theta point 

Can we identify one of the multicritical points discussed above with the usual theta or theta' p5|, 
f40| ] P| point? We recall that the theta point, which was solved in |Q involves a set of attractions 
between nearest and second nearest neighbours on the honeycomb lattice, with a partition function, 
for an even number of non contractible polymer chains, that was supposed so far to equal Z c [2/3\. 
The point k = 6 in our series has also g = 2/3, but a rather different operator content due to the Zq 
parafermions and the boundary conditions couplings. The exponent v = 4/7 is the same. For the 
theta point it is expected, and rather well established numerically, that the cross over exponent 
is tp = 3/7. This means that the perturbation physically identified as changing the attraction 
between monomers has dimension 

*-'-£-! ("> 

This is the dimension of the field X 5 after twisting, while for k = 6 the theory has an X s potential. 
Under this perturbation the theory flows to the k — 3 theory however, not the k — 1 theory! This 
can be interpreted in two ways. We can of course believe that the present set of multicritical points 
has nothing in common with the theta point, and that the possible concordance of some exponents 
is an accident. We can also try to minimize the available set of polymer systems by assuming 
that k = 6 is indeed the theta point, and that something was not understood in [^o). This would 

13 The distinction of theta and theta' points is a technical matter that we cannot discuss here. See Jio[ and 
subsequent comments and replies in Phys. Rev. Lett. 
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explain why th e numerical analysis of the problem is so difficult. Indeed in the present scheme, 
we can go from k — 6 to k — 1 with a crossover exponent ip = 3/7 at the k = 6 point only by 
passing through the k — 3 point, with v = 5/8. The mixture of these behaviours for the lattice 
polymers may indeed lead to very confusing measurments. We summarize some of the properties 
of the polymer multicritical series in figure 24. 



6 Conclusions 

In conclusion the two kinds of models we have introduced illustrate various aspects of N = 2 
theories. The geometrical reformulation of the vertex model leads to lattice analogs of the 
chiral ring. The multicritical polymers are not precisely defined yet But they provide a nice 
interpretation of the role of bosonic and fermionic degrees of freedom, a geometrical understanding 
of the boundary conditions, and hopefully physical realizations of N = 2 flow and spontaneous 
supersymmetry breaking. For both models, the bulk free energy can be readily obtained, without 
Bethe ansatz type computations. 

This work opens many physical questions. Among them, one would like to find an explicit 
N = 2 formulation of the Flory approximations, which maybe could be generalized to higher 
dimensions. Also, the series of multicritical polymer points seems a good place to apply the ideas 



of 1 45 . Notice that polymers can give rise to experiments in two dimensions |46| , and therefore it 
should be possible to observe in such systems consequences of N = 2 supersymmetry. 

On the more mathematical side, the study of interrelations between the singularity structure 
and the geometrical lattice models should be fruitful. 

Finally, it seems from our study that the most natural description of lattice models with N = 2 
continuum limit is geometrical. This is maybe not totally unexpected. For instance that the free 
energy can be obtained without computation, which is the kind of nice property one expects for a 
N = 2 model, seems to imply a geometrical setting where the non trivial observables are defined 
by topological properties, for instance the connectivity in polymers Maybe by pursuing this 
route one could indeed give lattice interpretations of the magic numbers that have been computed 
for topological theories |47], Q . 

Acknowledgments: We thank A.Leclair, P.Fendley, G.Moore, N.Read, N.Seiberg , N.Warner 
and A.Zamolodchikov for many useful discussions. 



We do not know their exact Boltzmann weights, as is usually the case for multicritical systems 
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Figure Captions 



Figure 1: The weights of vertices arc encoded in the R matrix. 
Figure 2: Graphical representation of the Temper ley Lieb algebra. 
Figure 3: The two vertices in the geometrical reformulation of the 6 vertex model. 
Figure 4: Geometrical representation of the symmetrizer on two strands. 
Figure 5: Geometrical representation of the fusion formula for R matrices. 
Figure 6: The various configurations that enter in the decomposition of the j = 1 fusion before 
applying the symmetrizers. 

Figure 7: The three remaining configurations after application of the symmetrizers. 
Figure 8: The value of the double contractible loop is [3] g . 

Figure 9: A typical fat graph entering the computation of the partition function. 

Figure 10: Factorization rule for one "ear" of the diagram. 

Figure 11: The same diagram after factorizing out the ears. 

Figure 12: Factorization of "handles". 

Figure 13: The four vertices of interaction for k = 3. 

Figure 14: Squares with e matrices on each face encode conveniently the algebraic form of the 
x l interaction terms. This figure corresponds to the interaction of figure 17. 
Figure 15: Diagonal propagation for the square lattice. 
Figure 16: Graphical computation of x 2 for k = 2. 
Figure 17: Formulation of x 3 in the Temper ley Lieb algebra. 
Figure 18: Strands propagating on a cylinder 

Figure 19: Insertion of x reduces the number of non contractible loops by two 
Figure 20: e matrix in the solid on solid basis. 

Figure 21: Configurations that contribute to the connected correlation function for the X 
operator. 

Figure 22: If the wiggly line represents anti periodic boundary conditions for fermions, a non 
contractible loop that crosses it once gets a factor 1 + 1 = 2, while a non contractible loop that 
crosses it twice gets 1 — 1 = 0. 

Figure 23: Multicritical polymers are obtained by adding multipieces interactions. 

Figure 24: Some features of the renormalization group flow in the multicritical polymers series. 
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